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MODIFICATIONS  TO  ITERATIVE  RECURSION  UNFOLDING  ALGORITHMS  AND 
COMPUTER  CODES  TO  FIND  MORE  APPROPRIATE  NEUTRON  SPECTRA 


INTRODUCTION 

The  unfolding  of  neutron  spectra  using  data  from  activation  foils,  Bonner 
spheres,  or  other  detectors  usually  involves  solving  an  equation  of  the  form 

E 

/■max 

Yj  =  I  A^E)  X(E)  dE  (1) 

Emin  j  =  1’2’3***M 

where  X(E)  represents  the  fluence  distribution  of  the  energy  of  the  neutrons, 

Y  is  the  response  of  the  jth  detector,  (E)  describes  the  response  of  the 
jth  detector  to  neutrons  of  energy  E,  and  M  is  the  total  number  of  detectors. 
Equation  (1)  is  a  degenerate  case  of  a  Fredholm  integral  equation  of  the  first 
kind. 

Although  several  methods  exist  for  the  formal  solution  of  first-order 
integral  equations,  none  of  these  methods  are  generally  applicable  when  the 
detector  response  function  A^ (E)  is  not  known  analytically.  This  is  the  case 
for  all  practical  systems  used  for  neutron  spectrometry.  In  practice,  A^ (E)  is 
experimentally  determined  and/or  calculated  and  is  usually  approximated  by  a 
response  matrix  having  discrete  values.  Equation  (.1)  is  then  replaced  by  a 
set  of  linear  equations,  and  one  is  left  with  the  problem  of  solving  M  equations 
in  N  unknowns,  where  M  is  the  number  of  detectors  available  and  N  is  the  number 
of  points  needed  to  define  the  neutron  spectrum.  In  matrix  notation  we  have 

Y  =  AX  .  (2) 

Usually  N  is  greater  than  M,  hence  no  unique  solution  exists,  and  even  when 
N  ^  M  the  response  matrix  A  is  usually  ill-conditioned  resulting  in  wildly 
oscillating,  sometimes  negative,  solutions  having  little  physical  significance. 
This  led  Gold  [1]  to  introduce  the  terms  "exact,",  "approximate",  and 
Manuscript  approved  March  2,  1984. 
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"appropriate"  to  characterize  the  solutions  to  equation  (2).  An  exact  solution, 
if  one  exists,  satisfies  the  equation  exactly  but  may  be  negative  or 
oscillatory.  Approximate  solutions  satisfy  the  equation  only  within  reasonable 
error  limits.  From  the  approximate  solutions,  selection  of  the  most  physically 
acceptable  solution  yields  an  appropriate  solution,  which  usually  is  not  unique. 

Methods  for  solving  equations  (1)  and  (2)  for  neutron  spectrometry  have 
been  summarized  by  Nachtigall  and  Burger  [2]  and  by  Patterson  and  Thomas  [3]. 

One  method  which  finds  a  non-negative  solution  to  equation  (2)  by  minimizing 
through  an  iterative  recursion  procedure  the  deviation  between  the  measured 
and  calculated  detector  responses  has  been  described  by  Scofield  [4]  and  Gold 
[1].  This  method  was  modified  by  O'Brien  et  al  [5]  and  Sanna  [6],  who  called 
the  computer  code  used  to  unfold  the  spectrum  B0N31G.  More  recently, 

Dorashenkc-  et  al  [7]  have  described  another  iterative  recursion  method  for 
spectrum  unfolding.  A  computer  code  SPUNIT,  which  uses  this  algorithm,  has 
been  written  at  Pacific  Northwest  Labs  by  Brackenbush  and  Scherpelz  [8] .  In 
addition  to  non-negativity  of  the  solution,  advantages  usually  listed  for  such 
iterative  recursion  methods  of  spectrum  unfolding  are  ease  of  programming,  and 
the  possibility  of  adding  smoothing  to  the  solution.  However,  the  application 
of  other  prior  knowledge,  such  as  cutoff  energy  or  preferred  spectral  shape, 
has  not  been  incorporated  into  the  solution.  It  is  our  purpose  to  show  how 
this  may  be  accomplished  in  order  to  find  more  appropriate  neutron  spectra 
with  these  solution  methods. 

EFFECT  OF  INITIAL  SOLUTION  AND  SMOOTHING  METHOD 

In  Fig.  1  we  show  a  typical  neutron  spectrum  that  has  been  unfolded  from 
Bonner  sphere  data  using  B0N31G  or  SPUNIT  and  the  Sanna  response  matrix  [9]. 

This  spectrum  was  produced  by  moderating  a  Cf-252  fission  source  by  placing 
it  at  the  center  of  a  60  cm  dia.  steel  sphere.  The  spectrum  was  further 
moderated  by  interposing  a  60  cm  x  60  cm  x  3.8  cm  thick  Lucite  slab  between 
the  steel  ball  and  the  Bonner  sphere  detectors.  Detector  counts  were  greater 
than  10K  giving  counting  errors  of  approximately  17,.  The  spectrum  shown  in 
Fig.  1  is  typical  of  spectra  unfolded  from  Bonner  sphere  data  in  that  there 
are  few  neutrons  in  the  intermediate  energy  region  (10-6  -  10-2  MeV)  giving 
rise  to  the  so-called  "Bonner  dip".  Also,  the  high  energy  end  of  the  spectrum 
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does  not  decrease  as  rapidly  as  calculations  suggest  [10].  Similarly  shaped 
spectra  are  often  also  obtained  with  other  unfolding  codes  such  as  LOUHI  [11] 
which  can  incorporate  prior  knowledge  into  the  solution  [12,  13].  We  shall 
see  that  the  physically  unreasonable  spectrum  shape  shown  in  Fig.  1  is  caused 
by  the  starting  solution  and  the  smoothing  algorithm  used. 

It  is  usually  stated  that  the  shape  of  the  initial  solution  has  little 
or  no  influence  on  the  final  solution  for  recursion  methods  [6,  7].  That 
this  is  not  the  case  is  illustrated  in  Fig.  2.  Using  SPUNIT  and  the  same  data 
as  in  Fig.  1,  spectra  were  unfolded  using  two  drastically  different  initial 
solutions.  One  initial  solution  was  a  1/E  spectrum,  i.e.,  a  straight  hori¬ 
zontal  line  on  the  lethargy  scale  shown.  The  second  had  alternate  initial 
spectral  values  differing  by  a  factor  of  10.  No  smoothing  was  used  in  either 
case.  Note  that  after  1000  iterations  the  spectrum  unfolded  with  alternate 
factor  of  10  initial  values  still  retains  this  characteristic.  Even  after 
5000  iterations  this  characteristic  remained.  Both  spectra  shown  in  Fig.  2 
fit  the  experimental  data  with  an  average  error  less  than  1%.  The  fact  that 
both  spectra  fit  the  data  so  well  indicates  that  Bonner  sphere  data  cannot 
produce  spectra  having  great  resolution,  hence  their  dependence  on  the  initial 
solution.  Using  B0N31G,  we  obtained  similar  results,  however,  3000  iterations 
were  necessary  to  fit  the  experimental  data  within  1%. 

The  smoothing  algorithm  also  affects  the  shape  of  the  spectrum  shown  in 
Fig.  1.  The  smoothing  used  was  as  follows: 

X  (i)  =  (SX(i-l)  +  X(i)  +  SX(i+l))/l+2S),  i  =  2,3*  “N-l  (3) 

s 

Xs(l)  =  X(l) ,  (3a) 

Xg(N)  =  X(N) ,  (3b) 

where  X  is  the  value  of  spectrum  before  smoothing,  Xg  is  the  value  of  the 
spectrum  after  smoothing,  and  S  is  the  smoothing  factor.  Smoothing  was  done 
on  each  iteration  with  S  =  0.01.  Sanna  [6]  used  an  equivalent  equation  with 
S  =  0.05,  smoothing  on  alternate  iterations.  X(l)  and  X(N)  were  also  smoothed. 

An  examination  of  equation  (3)  reveals  that  this  smoothing  also  tends  to 
bias  the  solution  to  a  1/E  spectrum.  To  eliminate  this  bias,  Sanna  later 
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removed  smoothing  completely  [14,  15].  Hence,  his  earlier  finding  [6]  that 
the  final  solution  was  independent  of  starting  solution  was  probably  caused 
by  the  introduction  of  smoothing.  Using  smoothing,  we  also  obtained  the  same 
final  solution  for  all  starting  spectra. 

BIASING  THE  SOLUTION  TO  OTHER  STARTING  SPECTRA 

The  previous  results  indicate  that  the  starting  solution  does  affect  the 
final  solution,  hence,  it  would  seem  appropriate  to  choose  the  most  physically 
reasonable  spectrum  for  a  starting  solution.  Most  neutron  spectra  encountered 
in  radiation  protection  work  are  produced  by  the  scattering,  moderation,  and 
absorption  of  neutrons  originally  produced  by  nuclear  fission,  particle  ac¬ 
celerators,  or  by  the  a,n  reaction  from  radioactive  sources.  These  processes 
tend  to  produce  spectra  that  can  be  characterized  as  having  a  high  energy  peak 
corresponding,  or  reduced  somewhat  by  moderation,  to  the  original  neutron 
source  energy,  a  (1/E)X  intermediate  energy  component  produced  by  elastic 
scattering,  and  a  thermal  peak  whose  magnitude  is  determined  by  the  atomic 
number  of  the  shielding  and  scattering  material  and  by  the  thermal  neutron 
absorption  cross-section  of  these  materials.  We  therefore  developed  an 
algorithm  called  MAXIET  [16]  for  use  in  YOGI  [17]  to  find  the  best  solution 
made  up  of  these  three  spectrum  components.  This  spectrum  can  then  be  used 
as  an  initial  solution.  Using  the  same  data  as  Fig.  1  and  2,  the  spectra 
obtained  using  only  MAXIET,  and  using  MAXIET  for  an  initial  solution  for 
SPUNIT,  are  shown  in  Fig.  3.  Using  only  MAXIET,  the  average  error  on  the 
detector  responses  is  1.1%;  using  MAXIET  plus  SPUNIT  this  error  is  reduced 
to  0.5%.  No  smoothing  was  used  in  either  case.  Note  that  using  MAXIET  and 
SPUNIT  gives  a  spectrum  that  fits  the  data  better,  and  agrees  better  with 
calculated  spectra  [10],  than  is  obtained  using  a  1/E  initial  solution. 

SMOOTHING  TO  OTHER  STARTING  SPECTRA 

As  previously  noted,  equation  (3)  biases  the  solution  to  a  straight  line, 
i.e.,  1/E  spectrum.  What  is  needed  is  some  method  to  bias  the  solution  toward 
any  chosen  spectrum,  e.g.,  the  MAXIET  initial  spectrum.  Perhaps  the  easiest 
way  to  do  this  is  to  make  a  transformation  on  the  response  matrix  and  the 
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initial  spectrum  so  that  the  transformed  initial  spectrum  Ls  a  straight  line, 
i.e. , 


AT(i,j)  =  A(i,j)  Xj (i)  i  =  1,2*  •  *N  (4) 

j  =  1,2* • *M 

X*(i)  =1  i  =  1,2* • *N  (5) 

T 

where  A  is  the  response  matrix,  A  is  the  transformed  response  matrix,  X  is 

T  L 

the  initial  spectrum,  and  Xj.  is  the  transformed  initial  spectrum. 

T  T 

The  solution  can  then  proceed  using  X  and  A  to  find  the  transformed 
T 

solution  X  ,  using  appropriate  smoothing,  e.g.,  equation  (3).  However,  since 
the  thermal  neutron  cross-sections  of  the  shielding  materials  can  tend  to 
cause  a  very  large  or  very  small  thermal  neutron  component,  we  do  not  smooth 
either  the  first  or  second  energy  intervals,  i.e.,  equation  (3)  is  replaced  by 

Xg(i)  =  (SXT(i-l)  +  XT(i)  +  SXT(i+l))/(l+2S),  (6) 

i  =  3,4* • *N-1 


Xg<l)  =  XT(1) 

(6a) 

XgC)  =  XT(2) 

(6b) 

Xg (N)  =  XT(N) 

(6c) 

Alternatively,  if  one  wishes  to  tie  the  Nth  spectral  value  to  a  fixed  (N+l)th 
value,  this  is  easily  accomplished  by  letting  the  index  in  equation  (6)  go  to  N 
and  eliminating  equation  (6c) . 

T 

After  solving  for  X  ,  the  final  solution  is  obtained  by  making  the  inverse 
transformation,  i.e.. 


x(i)  =  XI(i)XT(i),  i  =  1,2* • *N,  (7) 

To  illustrate  the  use  of  smoothing  to  different  starting  spectra,  errors  of 
approximately  5%  S.D.  were  introduced,  using  the  algorithm  in  [6],  to  the 
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data  used  in  the  previous  figures  to  simulate  larger  experimental  errors. 
Spectra  were  then  unfolded  using  a  1/E  and  a  MAXIET  starting  spectrum  using 
no  smoothing,  or  equation  (6)  for  smoothing  with  S  =  0.05.  These  results  are 
shown  in  Figs.  4  and  5.  Note  that  using  the  MAXIET  starting  spectrum  with 
smoothing  to  that  spectrum  gives  more  physically  reasonable  results  than  does 
the  1/E  starting  spectrum.  More  importantly,  it  can  be  seen  that  making  the 
matrix  transformation  allows  smoothing  to  any  chosen  initial  spectrum. 

PROGRAM  BUNKI 

The  modifications  which  we  have  described  to  find  more  appropriate 
neutron  spectra  using  B0N31G  or  SPUNIT  are  incoporated  in  a  computer  program 
called  BUNKI  which  is  listed  in  Appendix  A.  BUNKI  is  written  in  FORTRAN  IV 
and  was  coded  for  use  on  a  DEC-10  computer.  Tapes  of  the  program  and  the 
response  matrix  file  GIANT  can  be  generated  which  are  compatible  with  most 
mini  or  mainframe  computers.  Contact  the  authors  for  details. 

The  meanings  cfmost  of  the  variables  used  in  BUNKI  are  apparent  from 
their  use  in  the  program.  The  definitions  of  the  parameters  which  affect  the 
spectrum  and  which  may  be  input  by  the  user  are  listed  in  Appendix  B.  Those 
which  are  not  usually  changed  are  written  directly  into  the  program,  but  they 
could  also  be  input  by  the  user,  for  example,  at  lines  27400-28000.  Numbers 
in  parentheses  refer  to  relevant  line  number(s)  for  the  parameters. 

BUNKI  is  written  to  be  run  interactively  from  a  terminal.  Appendix  C 
is  an  example  of  a  BUNKI  session  in  which  two  spectra  are  generated  using 
the  same  Bonner  sphere  input  data.  In  the  first  case,  the  user  inputs  the 
initial  spectrum;  in  the  second,  the  MAXIET  algorithm  is  used  to  find  the 
initial  spectrum.  The  unformatted  variables  written  back  to  the  terminal 
as  a  guide  for  the  user  are  listed  at  lines  39200  and  59400  of  BUNKI.  The 
questions  which  BUNKI  asks  the  user  indicate  the  many  options  available  such 
as:  choice  of  detectors,  choice  of  unfolding  code,  choice  of  response  matrix, 
choice  of  initial  spectrum,  and  choice  of  the  number  of  energy  intervals,  in 
addition  to  the  input  parameters  in  Appendix  B. 

The  results  of  the  spectrum  unfolding  are  written  directly  to  disc  in 
a  file  called  FOR01.DAT  which  is  automatically  generated  by  the  DEC-10.  In 
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addition,  the  lethargy  spectrum  data  may  be  written  to  a  file  called  SPECX.DAT 
which  is  used  in  a  separate  program  to  create  plots  of  the  spectra.  The  output 
of  the  session  listed  in  Appendix  C  is  listed  in  Appendix  D. 

CONCLUSION 

By  properly  choosing  the  initial  solution,  and  by  smoothing  to  that 
solution  by  means  of  a  transformation  on  the  detector  response  matrix,  it  is 
possible  to  incorporate  prior  knowledge  into  the  neutron  spectra  obtained 
using  iterative  recursion  unfolding  algorithms.  We  have  written  a  computer 
code  called  BUNKI  which  incorporates  these  modifications  into  B0N31G  and 
SPUNIT.  These  modifications  should  also  be  useful  for  use  in  other  unfolding 
algorithms  and  computer  codes. 
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Neutron  spectra  of  Cf-252  moderated  by  30  cm  of  steel  and  3-81  cm 
of  Lucite.  The  spectra  were  unfolded  using  the  MAXIET  algorithm 
(□),  and  using  the  MAXIET  spectrum  as  an  initial  solution  for 
SPUNIT  (1000  iterations)  (+) .  No  smoothing  was  used  in  either 
case. 


REFERENCES 


1.  Gold,  R. ,  "An  Iterative  Unfolding  Method  for  Response  Matrices",  Argonne 
National  Laboratory,  Argonne,  IL,  ANL-6984  (1964). 

2.  Nachtigall,  D.  and  Burger,  G.,  "Dose  Equivalent  Determinations  in  Neutron 
Fields  by  Means  of  Moderator  Techniques",  in:  Topics  in  Radiation 
Dosimetry,  (Edited  by  F.  H.  Attix,  W.  C.  Roesch  and  E.  Tochlin)  pp .  385- 
459  (New  York:  Academic  Press)  (1972). 

3.  Patterson,  H.  W. ,  and  Thomas,  R.  H.,  Accelerator  Health  Physics  (New  York: 
Academic  Press)  (1973). 

4.  Scofield,  N.  E.,  "A  Technique  for  Unfolding  Gamma-Ray  Scintillation 
Spectrometer  Pulse  Height  Distributions",  U.  S.  Naval  Radiological 
Defense  Laboratory  San  Francisco,  CA,  USNRDL-TR-447  (1960). 

5.  O'Brien,  K. ;  Sanna,  R.;  and  McLaughlin,  J.;  "Inference  of  Accelerator 
Stray  Neutron  Spectra  from  Various  Measurements",  in:  First  Symposium 

on  Accelerator  Radition  Dosimetry  and  Experience,  USAEC-CONF-651109  (1965). 

6.  Sanna,  R.  S.,  "Modification  of  an  Iterative  Code  for  Unfolding  Neutron 
Spectra  from  Multisphere  Data",  Health  and  Safety  Laboratory  Energy 
Research  and  Development  Administration  New  York,  HASL-311  (1976) . 

7.  Doroshenko,  J.  J.;  Kraitor,  S.  N.;  Kuznetsova,  T.V.;  Kushnereva,  K.  K. ; 
and  Leonov,  E.  S.,:  "New  Methods  for  Measuring  Neutron  Spectra  with 
Energy  from  0.4  eV  to  10  MeV  by  Track  and  Activation  Detectors",  Nucl ■ 

Tech.  33,  296-304  (1977). 

8.  Brackenbush,  L.  W.  and  Scherpelz,  R.  I.,  "SPUNIT,  A  Computer  Code  for 
Multisphere  Unfolding",  in:  Computer  Applications  in  Health  Physics, 
Proceedings  of  the  Health  Physics  Topical  Meeting,  Pasco,  WA  (1984). 

9.  Sanna,  R.  S.,  "Thirty  One  Group  Response  Matrices  for  the  Multisphere 
Neutron  Spectrometer  over  the  Energy  Ranges  Thermal  to  400  MeV",  Health 
and  Safety  Laboratory,  U.  S.  Atomic  Energy  Commission,  NY,  HASL-767  (1973). 

10.  Ing,  H.  and  Makra,  S.,  Compendium  of  Neutron  Spectra  in  Criticality 
Accident  Dosimetry,  International  Atomic  Energy  Agency,  Vienna, 
STI/DOC/10/180  (1978). 

11.  Routti,  J.  T.  and  Sandberg,  J.  V. ,  "General  Purpose  Unfolding  Program 
LOUHI78  with  Linear  and  Nonlinear  Regularizations",  Comput.  Phys.  Commun. 

2A_,  119-144,  (1980). 

12.  Griffith,  R.  V.;  Slaughter,  D.  R.,;  Patterson,  H.  W.;  Beach,  J.  L.;  Frank, 

E.  G.;  Rueppel,  D.  W.,;  and  Fisher,  J.  C.;  "Multi-Technique  Characterization 
of  Neutron  Fields  from  Moderated  252Cf  and  239PuBe  Sources",  in: 

Proceedings  of  National  and  International  Standardization  of  Radiation 
Dosimetry .  International  Atomic  Energy  Agency,  Vienna,  IAEA  STI/PUB/471 
(1978). 


13 


13.  Endres,  G.W.R.;  Aldrich,  J.  M.;  Brackenbush,  L.  W.;  Faust,  L.  G.;  Griffith 
R.  V.;  and  Hankins,  D.  E.;  "Neutron  Dosimetry  at  Commercial  Nuclear  Plants", 
Pacific  Northwest  Laboratory,  Richland,  WA,  NUREG/CR-1769,  PNL-3585  (1981). 

14.  Sanna,  R.S.;  Hajnal,  F.;  McLaughlin,  J.E.;  Gulbin,  J.  F.;  and  Ryan,  R.  M.; 
"Neutron  Measurements  Inside  PWR  Containments",  Environmental  Measurements 
Laboratory,  U.  S.  Dept,  of  Energy,  New  York,  EML-379  (1980). 

15.  Sanna,  R.  S.,  Private  communication.  Environmental  Measurement  Laboratory, 
New  York  (1983). 

16.  Johnson,  T.  L.,  "A  Directed  Search  Algorithm  for  Finding  Physically 
Reasonable  Neutron  Spectra  from  Bonner  Sphere  Data"  to  be  presented  to 
the  Health  Physics  Society  Annual  Meeting,  New  Orleans,  LA,  June  1984. 

17.  Johnson,  T.  L.  and  Gorbics,  S.  G.,  "An  Iterative  Perturbation  Method  for 
Unfolding  Neutron  Spectra  from  Bonner  Sphere  Data",  Health  Phys.  41 , 

859  (1981). 


14 


APPENDIX  A 
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LISTING  OF  BUNKI  COMPUTER  CODE 


NOVEMBER  8,1983 


C  BUNKI 

C  THIS  PROGRAM  CALCULATES  THE  PARTICLE  FLUENCE,  DOSE, 

C  AND  DOSE  EQUIVALENT  SPECTRA  AS  A  FUNCTION  OF  NEUTRON 

C  ENERGY  USING  THE  B0N31G  OR  SPUNIT  UNFOLDING  CODES  TO 

C  SOLVE  THE  MATRIX  APPROXIMATION  OF  A  FREDHOLM  INTEGRAL 

C  EQUATION  OF  THE  FIRST  KIND.  THE  TOTAL  FLUENCE,  DOSE, 

C  DOSE  EQUIVALENT,  QUALITY  FACTOR,  AVERAGE  ENERGY,  AND 

C  RESPONSES  OF  SELECTED  DETECTORS  ARE  ALSO  CALCULATED. 

C  THE  i'nITIAL  SPECTRUM  MAY  BE  SPECIFIED  BY  THE  USER  OR 

C  A  1/E,  MAXWELLIAN,  INITIAL  SPECTRUM  MAY  BE  DETERMINED 

C  USING  MAXIET,  AN  ALGORITHM  ORIGINALLY  DEVISED  FOR  USE 

C  IN  YOGI,  THE  ITERATIVE  UNFOLDING  CODE  DEVELOPED  AT  THE 

C  NAVAL  RESEARCH  LABORATORY.  A  MATRIX  TRANSFORMATION 

C  ALLOWS  SMOOTHING  BIASED  TO  EITHER  INITIAL  SPECTRUM. 

C 

C  BUNKI  WAS  ORIGUNALLY  PROGRAMED  AT  THE  NAVAL  RESEARCH 

C  LABORATORY  IN  JULY,  1983,  BY  KIMBERLY  A.  LOWRY  AND 

C  TOMMY  L.  JOHNSON. 

C 

PROGRAM  BUNKI 
DOUBLE  PRECISION  BALL 

DIMENSION  ALETH( 1 2 , 3 1 ) , SPC( 31 ) , BCE( 1 2 ) , CTLD( 31), 

&  BCC( 12) ,CRAD(31) ,CREM(31 ) ,CNUTRK(31) ,HEAD( 20) , 

8  RAD( 31 ) , REM( 31 ) , EEND(32 ) ,CE( 3 1 ) , SPL( 32 ) , WDLETH( 31), 
&  CNTA(31 ) ,CA70(31 ) ,CHAN(31 ) ,SPLPLT(31 , 10) ,SPLI(32) , 

8  SPLMAX(32 ) ,LL(12) , BALL(12) , PCTERR( 12) ,CODE( 12) , 

5  PREM(31),ERRBCE(12),WHTBCE(12),A(108,31), 

6  SPLL(32) ,SS(32) ,BK(31 ,31 ) , VECT(31 ) 

C  ENERGY  INTERVAL  END  POINTS 

DATA  EEND/ 

S  1 .000E-08 , 4 . 140E-07 ,6 . 826E-07 , 1 . 445E-06 , 3 . 059E-06 , 

8  6 . 476E-06 ,1.371 E-05 , 2 . 902E-05 , 6 . 144E-05 , 1 . 301 E-04 , 

8  2 . 754E-04 , 5 . 929E-04 , I . 234E-03 ,2.6 1 3E-03 ,5.531 E-03 , 

8  1.1 71E-02 , 2 . 479E-02 , 5 . 247E-02 , 1 . 11  IE-01 ,2 . 237E-01 , 

8  4.508E-01 ,9. 072E-01 , 1 . 872E-00, 3 . 679E-00 , 7 . 408E-00 , 

8  1 .492E+01 , 2 . 58 1 E+01 , 4 .465E+01 , 7 . 725E+01 , 1 . 336E+02 , 

8  2.31 2E+02 , 4 . 000E+02/ 

C  NAMES  OF  THE  DETECTORS 
DATA  BALL/ 

8  8HBARE  ,8HBARE)-CD  ,SH2  INCH  ,8H2"+CD 

8  8H3  INCH  , 8H3"+CD  ,8(15  INCH  ,8H5"+CD 

8  8H8  INCH  ,8KI0  INCH  ,8H12  INCH  ,SH18  INCH  / 

C  FLUENCE  TO  DOSE  CONVERSION  FACTORS 
DATA  C RAD/ 

8  5 . 260E-10,6.085E-10,6 . 1 75E-10.6. 135E-10 , 6 .070E-10 , 

8  6.008E-10, 5.970E-10, 5.937E-10, 5.892E-10, 5 . 698E-10, 

8  5 .465E-10 , 5 . 251  F.-10 . 5. 149E-10, 5 .083E-10, 5 .039E-10 , 

8  5.629E-10.6.639E-10.7.347E-10, 1 .03SE-09, 1 .454F.-09, 

8  2 . 244E-09 , 3 . 406E-09 ,4 . 220F.-09 , 5 . 773E-09 , 6 . 622E-C9 , 

8  9 . 097E-09 , 9 . 610E-09 ,1.01 3E-08 , 1 . 1 37E-08 , 1 . 490E-08 , 

8  1 . 794E-08/ 

C  FLUENCE  TO  DOSE  EQUIVALENT  CONVERSION  FACTORS 


r 


05700 

DATA  CREM/ 

05800 

8 

1 . 1 51 E-09 , 1 . 232E-09 , 1 . 254E-09 , 1 . 244F.-09 , 1 . 227L-09 , 

05900 

8 

1 .211E-09, 1 . 195E-09, 1 . 179E-09, 1 . 161 E— 09 , 1 . 123E-09, 

06000 

& 

1 . 080E-09 . 1 . 040E-09 , 1 . 020E-09 . 1 . 008E-09 , 1 . 009E-09 , 

06100 

8 

1 . 514E-09 , 2 . 687E-09 . 4 . 7  75E-09 . 8 . 299E-09 ,1.4 1 6E-08 , 

06200 

8 

2 . 391 E-08 , 3 . 506E-08 . 3 . 995E-08 , 4 . 072E-08 . 4 . 11 1 E-08 , 

06300 

8 

4 . 356E-08 , 4 . 430E-08 .4.712E-03.5.026  E-08 , 5 . 5  36E-08 , 

06400 

& 

6.358E-08/ 

06500 

C 

FLUENCE  TO  NRL  TLD  READING  CONVERSION  FACTORS 

06600 

DATA  CTLD/ 

06700 

8 

06800 

8 

.4.. 4, .395, .37.. 33, .3. .27, .244, .22, .188, 

06900 

8 

.16,. 133, .11, .09, .07,. 05, 0,0, 0.0,0/ 

07000 

C 

FLUENCE  TO  'HANKINS'  TLD  READING  CONVERSION  FACTORS 

07100 

DATA  CHAN/ 

07200 

& 

.00007, 1.0, 1.0,. 9, .8, .75, .65, .6,. 5, .46, .42,.  3, 

07300 

& 

. 28,. 27,. 25,. 24,. 23,. 2,. 17,. 13,. 1,. 063,. 05,. 028,. 01 7,, 

07400 

8 

0,0, 0,0,0/ 

07500 

C 

DOSE  EQUIVALENT  TO  NEITRAK  144  READING  CONVERSION  FACTORS 

07600 

DATA  CNUTRK/ 

07700 

& 

0,0, 0,0. 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 

07800 

8 

.33,.  32.. 28, .35, .66,1.0, .57,0.2,0,0,0,0,0/ 

07900 

C 

DOSE  EQUIVALENT  TO  NTA  CONVERSION  FACTORS 

08000 

DATA  CNTA/ 

08100 

& 

0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0, 

08200 

8 

2.0,4.5,5.0,9.0,14.1,12.0,7.0,0,0,0,0,0/ 

08300 

C 

FLUENCE  TO  ANPDR-70  CONVERSION  FACTORS 

08400 

DATA  CA70/ 

08500 

8 

.0047 , .0057 , . 0058 ,.0061,. 0063 , . 0065 , . 0067 , .0070 , . 0072 , 

08600 

8 

.0075, .0077,. 0080, .0082, .0084 ,. 0081 ,. 0117, . 0169, 

08700 

8 

.0269,. 0386,. 0620,. 1152,. 1584, .2028, .1195. .0459, .0050, 

08800 

8 

0.0. 0,0.0/ 

08900 

C 

READ 

RESPONSE  MATRICES  FROM  FILE  'GIANT' 

09000 

OPEN  (UNIT-1,  FILE- 'GIANT') 

09100 

READ(1,*)A 

09200 

CLOSECUNIT-1 , FILE- 'GIANT' ) 

09300 

C 

COMPUTE  AVERAGE  ENERGY.  AND  LOG  WIDTH  OF  ENERGY  INTERVALS 

09400 

DO  1000  1-1,31 

09500 

C£(I)»( EEND( I)*EEND( 1+1 ))**.5 

09600 

1000 

WDLETH(I)-ALOG10(EEXD(I+l))-AL0C10(EEND(I)) 

09700 

c 

INPUT  FIXED  INITIAL  CONDITIONS  AND  DATA 

09800 

SLOPEJ-O 

09900 

PERSLP- .01 

10000 

THERM J-l .0 

10100 

THMMIN-0. 1 

10200 

THMMAX-10.0 

10300 

DEAD- . 000000 ! 32 

10400 

SHP-.01 

10500 

TSTRAT-.999 

10600 

c 

SET  OTHER  INITIAL  VALUES 

10700 

JX-1 

10800 

KX-1 

10900 

LX-1 

11000 

JJJ-0 

11100 

SPMX-1.0 

11200 

c 

INPUT  NUMBER  OF  ENERGY  INTERVALS 

11300 

1010 

WRITE(5, 1020) 

11400 

1020 

FORMAT/ ’  NUMBER  OF  ENERGY  INTERVALS?’) 

11500 

READ(5,*)JJ 

11600 

I F( J J . GE . J J J )CHNUM- ' Y ' 

11700 

JJJ-JJ+1 

11800 

MM-JJ 

11900 

C  INPUT  NUMBER  AND  KIND  OF  DETECTORS 

12000 

1030 

IF(JX.EQ.1)C0  TO  1060 

12100 

WRITE! 5 ,1040) 

12200 

1040 

FORMAT!'  CHANGE  DETECTORS?') 

12300 

1050 

FORMAT! Al) 

12400 

READ! 5 , 1050)CHDET 

12500 

IF!CHDET.EQ.'Y')G0  TO  1060 

12600 

IF(CHDET.EQ. 'N’)C0  TO  1140 

12700 

GO  TO  1030 

12800 

1060 

WRITE! 5, 1070) 

12900 

1070 

FORMAT!'  NUMBER  OF  DETECTORS?') 

13000 

REAC(5,*)KK 

13100 

1080 

WRITE! 5, 1090) 

13200 

1090 

FORMAT!’  type  DETECTOR  CODES!?  FOR  HELP)') 

13300 

READ! 5 , 1 100)C0DE( 1 ) 

13400 

1100 

FORMAT !A4) 

13500 

lf( CODEC 1).EQ. ’?')G0  TO  1120 

13600 

READ! 5 , 1 100 ) !CODE( I ) , I =2 , KK) 

13700 

DO  1110  I-l.KK 

13800 

LL(I)-0 

13900 

IF(C0DE(I).EQ.’0’)LL(I).1 

14000 

IF(C0DE(I).£Q.’0C')LL(I)-2 

14100 

IF(C0DE(I).EQ.'2’)LL(I)-3 

14200 

IF(C0DE(I).EQ.'2C')LL(I)«4 

14300 

IF(CODE(I).EQ.'3')LL(I)-5 

14400 

IF(C0DE(I).EQ. '3C')LL!I)=6 

14500 

IF(CODE(I).EQ.'5')LL(I)-7 

14600 

IF(C0DE(I).EQ.'5C')LL!I)-8 

14700 

IF(C0DE(I).EQ.'8')LL(I)=9 

14800 

IF(CODE(I).EQ. '10')LL(I)-10 

14900 

IF(C0DE(I ) .EQ . ' 12 ' ) LL( I )= 1 1 

15000 

IF(C0DE( I) .EQ. ' 18 ' )LL( t)=l  2 

15100 

IF(LL(I).EQ.O)GO  TO  1080 

15200 

1110 

CONTINUE 

15300 

GO  TO  1140 

15400 

1120 

WRITE(5, 1130) 

15500 

1130 

FORMAT!  lOX.’O . BARE’,/, 

15600 

8 

lOX.’OC . BARE+CADMIUM  COVER’,/, 

15700 

& 

10X,’2 . 2  INCH  BALL',/, 

15800 

8 

10X,’2C . 2  INCH  CADMIUM  COVERD  BALL',/, 

15900 

8 

10X/3 . 3  INCH  BALL',/, 

16000 

8 

10X,'3C . 3  INCH  CADMIUM  COVERED  BALL',/ 

16100 

8 

10X,  '5 . 5  INCH  BALL',/, 

16200 

8 

10X ,  ' 5C . 5  INCH  CADMIUM  COVERD  BALL’,/, 

16300 

8 

10X,  '8 . 8  INCH  BALL'  ,/, 

16400 

8 

lOX.'lO . 10  INCH  BALL',/, 

16500 

a 

10X , '  12 . 12  INCH  BALL',/, 

16600 

8 

10X, '  18 . 18  INCH  BALL’,/,/) 

16700 

GO  TO  1080 

16800 

C  SELECT  RESPONSE  MATRIX 

16900 

1140 

WRITE! 5,1145) 

17000 

1145 

FORMAT!'  TYPE  MATRIX  NAME!?  FOR  HELP)') 

17100 

READ( 5 , 1 100)RMTX 

17200 

IF(RMTX.EQ. '?' )G0  TO  1150 

17300 

LLL=1 

17400 

IF!RMTX.EQ.'SAN4')LLL-0 

17500 

IF!RMTX.EQ. ' SAN 1 ’ JLLL-12 

17600 

I F( RMTX . EQ . ' M60 ' ) LLL«24 

17700 

IF(RMTX.EQ. 'M65 ' )LLL*36 

17800 

IF(RMTX.EQ. 'MS13 ’ )LLL*48 

17900 

IF( RMTX . EQ . ' LOGN ' ) LLL-60 

18000 
18100 
18200 
18300 
18400 
18500 
18600 
18700 
18800 
18900 
19000 
19100 
19200 
19300 
19400 
19500 
19600 
19700 
19800 
19900 
20000 
20100 
20200 
20300 
20400 
20500 
206 00 
20700 
20800 
20900 
21000 
21100 
21200 
21300 
21400 
21500 
21600 
21700 
21800 
21900 
22000 
22100 
22200 
22300 
22400 
22500 
22600 
22700 
22800 
22900 
23000 
23100 
23200 
23300 
23400 
23500 
23600 
23700 
23800 
23900 


IF(RMTX.EQ. 'BARC')LLL=72 
IF(R>ITX.EQ.  ' UTA4 ' )LLL=84 
IF(  R>1TX .  EQ . '  UTA 1 ' )  LLL=96 
IF(LLL.EQ. 1 )G0  TO  1140 
GO  TO  1160 

1150  WRXTE(5, 1155) 

1155  FORMAT( 10X , ' SAN4 . SANNA  4MM  X  4MM  ' ,/, 

&  10X/SAN13 . SANNA  13MM  X  13MM  ',/, 

&  10X/M60 . M60 ' , / , 

&  10X,  'M65 . M65 ' ,  / , 

&  10X, 'MS 13 . MODIFIED  SANNA  13MMX13MM' ,/, 

5  lOX.'LOCNM . LOGNORMAL  '  J , 

6  lOX.'BARC . BHABA  MONTE  CARLO  ’/. 

&  10X,  'UTA4 . 4MM  X  4MM  U  OF  TEXAS  ',/, 

&  10X, 'UTA13 . 13MM  X  13MM  U  OF  TEXAS  ',/) 

GO  TO  1140 

C  SELECT  UNFOLDING  CODE 
1160  Will  TE(5, 1165)  ■ 

1165  FORMATC  TYPE  UNFOLDING  C0DE(?  FOR  HELP)’) 

READ( 5 , 1 1 00 ) UN  FOLD 

IF( UNFOLD. EQ.’?')GO  TO  1170 

IF(UNF0LD.EQ.'B0N3'. OR. UNFOLD. EQ. ’SPUN’ )G0  TO  1180 
GO  TO  1160 

1170  WRITE(5, 1175) 

1175  FORMAT! 10X, 'B0N31G . HASL  UNFOLDING  CODE’./, 

&  10X . ' SPUNIT . PNL-SOVIET  UNFOLDING  CODE '  ,  /  ) 

GO  TO  1160 
1180  CONTINUE 

C  SELECT  APPROPRIATE  COLUMNS  FROM  MATRIX  AND  ADJUST 
C  MATRIX  TO  UNFOLD  PER  UNIT  LETHARGY 
DO  1190  1=1, KK 

L=LL(I)+LLL 
DO  1190  J=1,JJJ 

1190  ALETH(I.J)=A(L,J)*WDLETH(J) 

C  INPUT  VARIABLE  INITIAL  CONDITIONS  AND  DATA 
1200  WRITE(5, 1210) 

1210  FORMAT! '  PUT  SPECTRA  IN  FILE  FOR  PLOTTING?’) 

READ! 5, 1050) FILE 

1220  IF(FILE.EQ. 'Y')GO  TO  1230 

IF(FILE.EQ.’N’)GO  TO  1250 
GO  TO  1200 

1230  WRITE(5, 1240) 

1240  FORMAT! '  FILE  NUMBER  OF  THE  SPECTRUM?  (1  TO  10)') 
READ(5,*)KX 

IF(KX.LT. 1 .0R.KX.GT.10)G0  TO  1230 
1250  IF(JX.EQ.1)G0  TO  1270 

WRITE! 5, 1260) 

1260  FORMATC  CHANGE  THE  HEADING?') 

READ! 5 , 1050)CHHEAD 
IF(CHHEAD.EQ.'Y')GO  TO  1270 
IF(CHHEAD.EQ.  '.V')CO  TO  1320 
GO  TO  1250 

1270  WRITE  (5,1280) 

1280  FORMATC  TYFE  THE  HEADING  (<80  CHARACTERS)') 

READ! 3 , 1290) (HEAD! I ) , 1=1 , 20) 

1290  FORMAT! 20A4) 

1320  WRITE! 5, 1330) 

1330  FORMATC  SEARCH  FOR  MAXWELLIAN,  1/E  INITIAL  SPECTRUM? ' ) 
READ! 5 , 10501SRHMA 

IF( SRHMA.SE. 'Y' .AND.SRHMA.NE. ’N')GO  TO  1320 
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1 


24000 

I F ( SRHMA . N  E . SRHMAX )  LX= 1 

24100 

SRHMAX-SRHMA 

24200 

IF( SRHMAX . EQ . ’ Y ' )G0  TO  1380 

24300 

IF(SRHMAX.EQ. 'N')CO  TO  1340 

24400 

GO  TO  1320 

24500 

1340 

IF(LX. EQ. 1 )G0  TO  1360 

24600 

IF(CHNUM.EQ. ' Y' )G0  TO  1360 

24700 

WRITE( 5, 1350) 

24800 

1350 

FORMAT ( '  CHANGE  INITIAL  SPECTRUM? ' ) 

24900 

READ( 5 , 1050  JCHISPC 

25000 

IF(CHISPC.EQ. ' Y ' )G0  TO  1360 

25100 

IF(CHISPC.EQ.'N’)GO  TO  1480 

25200 

GO  TO  1340 

25300 

1360 

W’RITE(5, 1370)JJJ 

25400 

1370 

FORMATC  TYPE  INITIAL  SPECTRUM  (',12,'  VALUES)') 

25500 

READ( 5 ,*)(SPLI( I) ,1=1 , J3J) 

25600 

GO  TO  14 80 

25700 

1380 

IF(LX . EQ. 1 )G0  TO  1400 

25800 

WRITE( 5 ,1390) 

25900 

1390 

FORMATC  CHANGE  MAXWELLIAN  TEMP,  SHAPE,  OR  PERTURBATION 

26000 

READ(5, 1050)CHTEMP 

26100 

IF(CHTEV1P.EQ.'Y')G0  TO  1400 

26200 

IF(CHTEMP.EQ. 'N')G0  TO  1480 

26300 

GO  TO  1380 

26400 

1400 

WRITE( 5 , 1410) 

26500 

1410 

FORMATC  TYPE  .MAXWELLIAN  TEMP,  SHAPE,  AND  PERTURBATION’; 

26600 

READ( 5 ,*)TEMPIJ , SHAPE, PERTMP 

26700 

1480 

IF(JX.EQ.1)G0  TO  1500 

26800 

WRITE( 5 , 1490) 

26900 

1490 

FORMATC  CHANGE  FIT  PARAMETERS?') 

27000 

READ( 5 , 1 050)CHFIT 

27100 

IF(CHFIT.EQ. ' Y ' )G0  TO  1500 

27200 

IF(CHFIT.EQ. 'N')GO  TO  1520 

27300 

GO  TO  1480 

27400 

1500 

WRITE(5, 1510) 

27500 

1510 

FORMAT ( '  TYPE:  END  TEST(T) , ' , 

27600 

& 

/,'  SMOOTHING  FACTOR,  ', 

27700 

& 

/,'  CALIBRATION  FACTOR,  ’, 

27800 

& 

/,’  ITERATIONS  BEFORE  ERROR  TEST,  ', 

27900 

a 

/, '  AND  MAXIMUM  NUMBER  OF  ITERATIONS') 

28000 

READ(5,  *  )TSTPER ,  S.MO ,  CAL ,  ITRTST ,  ITRMAX 

28100 

PERTKM=1.0+20*PERSLP 

28200 

PERE-1 .0+10*PERSLP 

28300 

IF( SRHMAX. EQ. 'N'. AND. ITRMAX. EQ.O)GO  TO  1820 

28400 

TSTPER=(  KK*TSTPER*-*2 )  / 10000 

28500 

1520 

TEMPI=TEMPIJ 

28600 

SLOPEI -SLOPE J 

28700 

THERMI-THERMJ 

28800 

IF(JX.EQ. 1 )G0  TO  1540 

28900 

IF(CHDET.EQ. ' Y ’ )G0  TO  1540 

29000 

1525 

WRITE(5, 1530) 

29100 

1530 

FORMATC  CHANGE  BALL  DATA?') 

29200 

READ(  5 , 1 050)CIIBCE 

29300 

IF(CHBCF..EQ.'Y')GO  TO  1540 

29400 

IF(CHBCE. EQ. ' N ' )C0  TO  1590 

29500 

CO  TO  1525 

29600 

1540 

DO  1560  1=1, KK 

29700 

L=LL( I ) 

29800 

WRITE(5, 1550)BALL(L) 

29900 

1550 

FORMATC  TYPE  ', AS,  ’BONNER  SPHERE  COUNT,  2  ERROR') 
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30000 

1560 

READ!  5 ,*) ( BCEC 1 ) .ERRBCFJ 1 1 ) 

30100 

C  SUM  ERRORS,  MAKE  DEAD  TIME  CORRECTION  TO  BALL  CPVNTS 

30200 

SUMWHT-0 

30300 

DO  1570  I-l.KK 

30400 

SU>!WHT-SUMVHT+ERRBCE(  I ) 

30500 

1570 

BCE( I )-ECE( I )/ ( 1 .0-BCE( I )*DEAD) 

3060  0 

C  CALCULATE  BALL  COUNT  ERROR  WEIGHTS 

30700 

DO  1580  I-l.KK 

30800 

1580 

WHTBCE( I ) -SUMWHT/ ( KK»ERRBCE( I ) ) 

30900 

C  CALCULATE  MAXWELLIAN,  1/E  SPECTRUM  IE  REQUIRED 

31000 

1590 

IF(SSHMAX.EQ.  '.N')CO  TO  1820 

31100 

C  BEGINNING  OF  MAXIET  ALGORITHM 

31200 

C 

INITIALIZE  FIT  PARAMETERS 

31300 

ERRORE- 123456789 

31400 

TEMP-TEMPI 

31500 

ERRORM-ERRORS 

31600 

1600 

SLOPE-SLOP El 

317CO 

THERM-THEKMI 

31800 

C 

CALCULATE  MAXWELLIAN  SPECTRUM 

31900 

SPMX-0 

32000 

DO  1610  I-l.JJJ 

32100 

SPLMAX(  I  )=(  CE(  I  )**  1 . 5  )*(  EXP(  -CE(  I  )/TEMP) ) 

32200 

IF(SPLMAX( I ) .GT.SPMX )SPMX-SPLMAX( I ) 

32300 

IF(SPMX .  EQ .  SPL>IAX(  1 )  )G0  TO  1610 

32400 

SPLMAX  ( 1  )-SP>V‘:**SHAPE*SPLHAX(  I )**(1 . 0-SHAPE) 

32500 

IF(SFLMAXd)  .LT.SPLMAX(I-1)*3HP)SPLMAX(I)=SPLMAX( 

32600 

1610 

CONTINUE 

32700 

C 

calculate  i/e  initial  spectrum 

32800 

ERROR-123456739 

32900 

ERRORE- 123456789 

33000 

HGTE=SPMX*PER£ 

33100 

1620 

ERRORT-ERROR 

3320Q 

1630 

HGTE-HGTE/PERE 

33300 

1640 

DO  1650  T-I.JJJ 

33400 

1650 

SPLI(  I  )=HGTE*CE(  I  )»»SLOPE 

33500 

C 

COMBINE  MAXWELLIAN  AND  1/E  SPECTRA 

33600 

DO  1660  1=1, JJJ 

33700 

IF(SPLI(I).LT.SPLMAX(I))CO  TO  1670 

33800 

SPLI(I)-(SPLI( I)+SPLMAX(I) )*0,5 

33900 

1660 

CONTINUE 

34000 

GO  TO  1630 

34100 

1670 

DO  1680  J-I.JJJ 

34200 

1680 

SPLI(J)-SPLMAX(J) 

34300 

C 

ADJUST  THERMAL  ENERGY  BIN 

34400 

SPLI ( 1 )-3PLI ( 2 )*THERM 

34500 

C 

CALCULATE  SPHERE  RESPONSES  AND  SUM  FROM  SPECTROt 

34600 

DO  1690  M«I ,.\K 

34700 

BCC( M)=0 

34800 

DO  1690  J-l.JJJ 

34900 

1690 

BCCI  M)-BCC(M  )+ALETll(  M .  J  )"'SPLT  ( J  ) 

35000 

C 

CALCULATE  SUMS  OF  SPHERE  DATA 

35100 

SUMBCE-0 

352CO 

SUMBCC-0 

35300 

DO  1700  I-l.KK 

35400 

SUMBCF.=SUMBCE+BCE(  I ) 

35300 

1700 

SUMBCC=S12'IECC+BCC(  1 ) 

35600 

C 

NORMALIZE  CALCULATED  SPHERE  RESPONSES 

35700 

C 

TO  EXPERIMENTAL  DATA 

35800 

RNOEM-SUMBCE/SUMBCC 

35900 

DO  1710  I-l.KK 

20 


36000 

36100 

36200 

363CO 

36400 

36500 

36600 

36700 

36800 

36900 

37000 

37100 

372GO 

37300 

37400 

37500 

37600 

37700 

37800 

37900 

38000 

38100 

38200 

38300 

38400 

38500 

38600 

38700 

38SC0 

38900 

39000 

39100 

39200 

39300 

39400 

39500 

39600 

39700 

39800 

39900 

40000 

40100 

40200 

40300 

40400 

40500 

40600 

40700 

40800 

40900 

41000 

41100 

41200 

41300 

41400 

41500 

41600 

41700 

41800 

41900 


1710  BCC ( I )=BCC(  I  )*RNORM 

C  CALCULATE  ERROR  ON'  FIT 
ERR0R=0 
DO  1720  1=1, KK 
ERR=( BCC( I )-BCE< I ) )/BCE( I ) 

1 720  ERROR=ERROR+WHTBCE( I)*ERR*ERR 

IF( ERROR .LT.ERRORT)GO  TO  1620 
HGTE=HGTE*PERE 

1730  IF(ERRORT.GE.ERRORE)GO  TO  1740 

C  SAVE  BEST  VALUES  OF  FIT  PARAMETERS 
ERRORE=ERRORT 
HGTE£=HCT£ 

THERME-THERM 

SLOPEE-SLOPE 

MX-0 

c  change  slope 

1 740  SLOPE=SLCPE+P£RSLP 

IF(MX.EQ.1)SL0PE«SL0PEE 
C  CHANGE  THERMAL  BIN 

TH  ERM=TH  ERME*  P  ERTKM 
I F( THERM . GE . THMMAX ) GO  TO  1750 
I F( MX . E Q . 0)THERM=THERME 
MX-MX+1 

IF(MX -GT. 10)G0  TO  1750 
HGTE=HGTE*P£RE* ( 1 . 0+ 1 0*PERSLP ) 

IF(MX . EQ. 1 )HGTE=HGTEE*PERE*PERTHM 
C  RESET  ERROR.  SEARCH  FOR  BETTER  FIT  PARAMETERS 

ERR0RT=1 23456789 
GO  TO  1640 

C  CALCULATE  ERROR  ON  FIT 
1750  PERR0R=100il(ERR0RE/KiC)**.5 
C  WRITE  BEST  VALUES  OF  FIT  PARAMETERS  TO  TERMINAL 
WRITE(5, »)TEMP, SHAPE, HGTEE, SLOPES, THERME, PERROR 
IF( ERRORS. GE.ERRORN)CO  TO  1760 
C  SAVE  BEST  VALUES  OF  FIT  PARAMETERS 

TEMPM=TEMP 
HGTEM=HGTEE 
THERMM-THERME 
SLOP£M*SLOPE£ 

ERRORM=ERRORE 

C  CHANGE  MAXWELLIAN  TEMPERATURE  IF  REQUIRED 
IF( PERTMP. F.Q .0)G0  TO  1760 
TEMP=TKMP-PERTMP 

C  RETURN  AND  SEARCH  FOR  BETTER  PARAMETERS  IF 
C  MAXWELLIAN  TEMP  IS  IN  RANGE 

IF(TEMP.CT.PERTMP)GO  TO  1600 
1760  CONTINUE 

C  IF  FINAL  PARAMETERS  EQUAL  INITIAL  PARAMETERS  CAH.VGE 
C  INITIAL  PARAMETERS  AND  CONTINUE  SEARCH 
A  1=0 

IF'  ( SLOPEM .  EQ .  SLOPEI )  A 1  *  1 
A2=0 

IF(THER'!M.  EQ.THERMI .  AND. THERMI .  GE.THMMIN )A2=1 
A3=0 

IFCTEMPM.F.O.  TEMPI.  AND.  TEMPM.LT.TEMPIJ+10*PERTMP)A3=1 

I F( PERTMP. EQ.0)A3=0 

IF(A1  .F.Q.  1  )SL0PEI=SLOPEI-10*PERSLP 

IF(A2 . EQ. 1 )THEKMI=THERMI/PESTHM**3 

IF(A3.EQ.  1  )TEMPI»TEMPI+3*PERTMP 

IF(A1+A2+A3.GT. .5)G0  TO  1590 
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42000 

C 

IF  NOT,  CALCULATE  INITIAL  SPECTRUM  WITH  BEST  PARAMETERS 

42100 

SPMX-0 

42200 

DO  1770  I-l.JJJ 

42300 

SPLMAX(  I  )=(CE(  I  )**1 . 5)*(EXP(-CE(  D/TEMPM) ) 

42400 

IF(SPLMAX( I ) .GT.SPMX )SPMX=SPLMAX( I ) 

42500 

IF( SPMX . EQ . SPLMAX( I ) )G0  TO  1770 

42600 

SPLMAX( I )=SPMX**SHAPE*SPLMAX( I Js*( 1 .0-SHAPE) 

42700 

IF(SPLMAX(  I)  ,LT.SPLMAX(I-1  )*SHP)SPLMAX( I)=SPLMAX(  1-1 )' 

42800 

1770 

CONTINUE 

42900 

DO  1780  I-l.JJJ 

43000 

1730 

SPLI ( I )-HGTEM*CE( I )**SLGPEM 

43100 

DO  1790  I-l.JJJ 

43200 

IF(SPLI(I).LT.SPLMAX(I))GO  TO  1800 

43300 

SPLI(I)=(SPLI(I)+SPLMAX(I))*O.S 

43400 

1790 

CONTINUE 

43500 

1800 

DO  1810  J-I.JJJ 

43600 

1810 

SPLI ( J ) -SPLMAX ( J ) 

43700 

SPLI ( 1 ) -SPLI ( 2 ) *THERMM 

43800 

C  COMPLETION  OF  MAXIET  ALGORITHM 

43900 

C  TRANSFORM  MATRIX  TO  CONSTANT  INITIAL  SPECTRUM 

44000 

1820 

DO  1825  I-I.JJJ 

44100 

DO  1825  K-l.KK 

44200 

1825 

ALETH( K , I )=ALETH(K . I )*SPLI ( I ) 

44300 

DO  1830  I-l.JJJ 

44400 

1830 

SPL(I)«1 

44500 

C  CALCULATE  SPHERE  RESPONSES  AND  SUM  FROM  INITIAL  SPECTRUM 

44600 

DO  1840  M-l.KK 

44700 

BCC(M)=0 

44800 

DO  1840  J-l.JJ 

44900 

1340 

BCC (M)=3CC(M) + ALETH (M,J)*SPL(J) 

45000 

IF(SRHMAX.EQ. ' N ' .AND.ITRMAX.EQ.O)GO  TO  2000 

45100 

C  NORMALIZE  CALCULATED  SPHERE  RESPONSES  AND  INITIAL 

45200 

C  SPECTRUM  TO  EXPERIMENTAL  DATA 

45300 

SUMECC-0 

45400 

SLMBCE-0 

45500 

DO  1850  1=1, KK 

45600 

SUMBCE-SUMBCE+BCEd) 

45700 

1850 

SUMBCC=SUMBCC+BCC( I ) 

45800 

RNOP.M-SUMECE/  SLMBCC 

45900 

DO  1860  1=1, KK 

46000 

1860 

BCC(I)=BCC(I)*R.VORM 

46100 

SUMBCC-SIM3CE 

46200 

DO  1870  I-l.JJJ 

46300 

SPL(  I)=SPL(  1} -R.VCRM 

46400 

1370 

CONTINUE 

46500 

C  CALCULATE  ERROR  ON  FIT 

46600 

ERROR-O 

46700 

DO  1880  1=1, KK 

46800 

ERR=(BCC(I)-BCE(I) )/BCE(I) 

46900 

1880 

ERROR-ERROR+WfTBCE( I ) *ERR*ERR 

47000 

1F( ITRMAX . EQ.O)GO  TO  2C00 

47100 

ITER-0 

47200 

1890 

ERRORU-ERROK 

47300 

C  SELECT  UNFOLDING  ALGORITHM 

47400 

IF( UNFOLD . EQ . ' E0N3 ' )G0  TO  1945 

47500 

C  BEGINNING  OF  SPUNIT  UNFOLDING  ALGORITHM 

47600 

IF( ITER ,GT.O)GO  TO  1905 

4  77CO 

DO  1900  J-l.JJ 

47800 

SS(J)=0 

47900 

DO  1900  I-l.KK 

4SU0U 

!«00 

SSf  J)-SS(  J>ALF.T!i(  I,J).'3CEl  I  ) 

48  ICO 

lr*05 

DO  194.)  K=1,ITKT-T 

48200 

ITER-lTF.R+1 

4830U 

DO  1910  J- 1 , JJ 

48400 

SFLL( J )=0 

48  5oO 

DO  1910  I-l.KK 

48600 

1910 

SP!.L(J)-SPEL(j)»(SPL(J)*ALETH(I,  J)/(S3(J;’  3CC{ 

48700 

DO  1920  J-3.MM 

48800 

1920 

SPL( J)-(SPL1  '  i  ‘  »SMO+SPLL( J)+SPLL( J+l )*SMO>/( 

48900 

SPLd)-SPLL(l) 

49000 

SPL(2)=SPLL(2) 

49100 

IF(MM.EQ. JJ)GO  TO  1925 

49200 

SPL( JJ)=SPLL( JJ) 

49300 

1925 

DO  1930  M=1,KK 

49400 

BCC(M)=0 

49500 

DO  1930  J.l.JJ 

49600 

1930 

BCC(M)-3CC(M)+ALETK( M, J )*SFL( 1 ) 

49700 

1940 

CONTINUE 

49800 

C  END 

OF"  SPUNIT  UNFOLDING  ALGCRIT1LM 

49900 

GO  TO  1990 

50000 

C  BEGINNING  OF  P0N31G  UNFOLDING  ALGORITHM 

50100 

1945 

IF( ITER.GT.OJGO  TO  I960 

50200 

DO  1950  1=1, JJ 

50300 

DO  1950  J=1 ,JJ 

50400 

BK ( J , I )  =0 

50500 

DO  1950  M= 1 . KK 

50600 

1950 

BK( J ,  I )=EK( J , I )+ALETH(M, J)*ALETH(M. I ) 

50700 

DO  1955  1=1, JJ 

5O800 

VECT( I )=0 

50900 

DO  1955  J»!,KK 

51000 

1955 

VECT( I )=VECTCI)tALETH( J , I )*BCE( J) 

51100 

1960 

DO  1930  N=i , ITRTST 

51200 

ITER=ITER+1 

51300 

DO  1970  J=1,JJ 

51400 

AX=0. 

51500 

DO  1965  M»1,JJ 

51600 

1965 

AX=SPL(M)*BK ( J , M )+AX 

51700 

1970 

SPLL( J)=SPL( J )*VECT( J )/ AX 

51800 

DO  1975  J=3,MM 

51900 

1975 

SPL(J)=(SPLL(J-l)»SMO+SPLL(J)+SPLL(J+l)*SMO)/( 

52000 

SPL( 1 )=SPLL( 1 ) 

52100 

SPL(2)=SPLL(2) 

52200 

IF(MM,EQ. JJ)GO  TO  1980 

52300 

SPL( JJ)=SPLL(JJ) 

52400 

19S0 

CONTINUE 

52500 

DO  1985  M=1  ,  KK 

52600 

BCC(M)=0 

52700 

DO  1985  J-l.JJ 

52800 

1985 

BCC(M)=BCC(M)+ALETH(M, J)*SPL(  J ) 

52900 

C  F.ND 

OF  B0N3IG  UNFOLDING  ALGORITHM 

53000 

1990 

CONTINUE 

53100 

C  CALCULATE  ERROR  ON  FIT 

53200 

ERROR =0 

53300 

DO  1995  1=1, KK 

53400 

ERR=( BCC( I )-BCE( I ))/BCE( I ) 

53500 

1995 

ERROR=EREOR+WHTBCE(  I  )  ,ERR',ERR 

53600 

C  TEST 

■  FOR  COMPLETION,  CONTINUE  IF  NOT  COMPLETE 

53700 

IF(  ERROR/' ERRORU.LT.TSTRAT.  AND.  ITER+ITRTST.LE.I' 

53300 

a 

•  AND. ERROR .GT.TSTPF.R ) GO  TO  1890 

53900 

C  IF  COMPLETE,  DO  INVERSE  TRANSFORM  OF  SPECTRUM  AND  M. 

>)) 

d«SMO) 


L2,fSMO) 


KMAX 

IRIX 


A 
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54000 
54100 
54200 
54300 
54400 
54  500 
546C0 
54700 
54800 
54900 
55000 
55100 
55200 
55300 
55400 
55500 
55600 
55700 
55800 
55900 
56000 
56100 
56200 
56300 
56400 
56500 
56600 
56700 
56300 
56900 
57000 
57100 
57200 
57300 
57400 
57500 
57600 
57700 
5  7  S00 
57900 
58000 
58100 
58200 
53300 
53400 
5S500 
5S600 
53700 
58800 
58900 
59000 
59100 
59200 
59200 
59400 
59500 
59600 
59700 
598G0 
59900 
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2000  DO  2005  I-l.JJJ 

200c‘  £TL(  I  )=SPL(  i ) 'SFLI( I ' 

DO  2010  1=1, KK 
DO  2010  J-l.JJJ 

2010  ALETH(1.J)-ALETH(I.J)/SPLI(J) 

c  calculate  output  values 

IF(SRHMAX.EQ. ’N* . AND. ITRMAX . EQ.OJGO  TO  2030 
2015  HGTEM-0 . 5*HCTEM/SPMX 

SUMERR-0 
DO  2020  1=1, KK 

PCTERRt I )-10C*( BCC( I )-BCE( I ) )/BCE(  I ) 

2020  SUMERR-Sl>tERR+PCTERR(  I )*PCTERR(  I ) 

PERR0R=(SL'MERR/KK)**.5 
2030  SUMSPC=0 

SUMRAD-0 
SUMREM-0 
SUMEXS-0 
SUMTLD-0 

sCmhan-o 

SUMNTR-0 

SUMNTA-0 

SUMA70-0 

DO  2035  1-1,31 

IF(I.GT.JJ)SPL(I)-0 

SPL(I)=SPL(I)*CAL 

SPLPLT ( I , KX ) -SPL( I ) 

SPC( I )-SPL( I )*WDL£TH( I ) 

SUMSPC-SUMSPC+SPC( I ) 

REM( I )«CREM( I )*SPC( I ) 

SUMREM-SUMREM+REX(1) 

RAD(I)«CRAD( I )*SPC( I ) 

SU>!RAD=SUMRAD+RAD(  I ) 

SUMEXS-SUMF.XS+CE  ( I  )*SPC(  I ) 

SU>mD-SUMTLD+CTLD(  I )  *SPC(  I ) 

SUMHAN-SUMHAN+CHAN ( I )*SPC( I ) 

SUMNTR=SUMNTR+CNUTRK( I )*REM( I ) 
SUM.VTA=SLMA'TA+CVT.4(I)*REMf  I) 

2035  SUMA70-SUMA  70+CA  70( I )*SPC( I ) 

QF-SUMREM/SUMRAD 

AVEEN-(SU>!EXS-CE(  1  )*SPC(  1 )  )/(SUMSPC-SPC(  1 )) 

SUMTLD- ( SLMTLD/SUMREM  )/4 . 064E+06 
SUMHAN-(  SUMHA.V/ SUMREN )  /  2 . 085E+06 
SUMNTR«( SUMNTR/SUMR EM ) / . 5748 
SUMNTA-( SUMNTA/SUMREM ) /8 . 033 
SUMA70=(SUMA70/SUMR£M)/4 . 755E+06 
DO  2040  1=1, JJ 

2040  PREM(  I  )=100*(REM(  D/SUMREM) 

C  SET  UNCOMPUTED  PARAMETERS  TO  0,  WRITE  PARAMETERS  TO  TERMINAL 
IF(SRHMAX.EQ. ’ N ’ )THERMM=0 
I F( SRHMAX . EQ . ’ N ’ )TEMPM-0 
IF(SRIIMAX.EQ.  fN"  )HGTEM«0 
IF(SRHMAX.EQ. 'N ' )SLOPEM=0 
IF( SRHMAX. EQ.  ’?!'  .AND.ITRMAX.EQ.O)PERROR=Q 
IF( ITRMAX . EQ.O) ITER=Q 

WRITE( 5 , *)TEMPM , SHAPE . HG7EM , SLOPEM . THERMM , PERROR .ITER 
C  DECIDE  IF  RESULTS  ARE  WORTH  KEEPING 
2043  WRITE( 5, 2050) 

2050  FORMATC  SAVE  THESE  RESULTS?') 

READ( 5 , 1030)SAVE 
IF(SAVE.EQ. 'N')GO  TO  2190 
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60000  IF(SAVE.EQ.  'Y')G0  TO  2055 

60100  CO  TO  2045 

60200  C  WRITE  RESULTS  TO  DISC 
60300  2055  WRITE! 1 ,2060)(HEAD(I) ,I»l ,20) 

60400  2060  FORMAT! ' 1 ' , 20A4 ) 

60500  IF( ITRMAX . EQ.O)CO  TO  2030 

60600  WRITE( 1 , 2070 ) RMTX . USFOLD , TEMPM . SHAPE , CAL , SMO , PERROR , ITER 

60700  2070  FORMAT! IX, / ,4X, 'RESPONSE  UNFOLD  MAXWELL  CALIB. 

60800  &  SMOOTH  PERCENT  NO.  OF',/. 

60900  i  4X, 'MATRIX  CODE  TEMP, SHAPE  FACTOR  FACTOR  ERROR 

61000  a  ITERATIONS', 

61100  S  / , IX ,A3, A10.FS.2 , ' , ' ,F4.2,F8.4,2! F10.4) , IS) 

61200  GO  TO  2100 

61300  2080  WRITE! 1 ,2090)  RMTX , TEMPM, SHAPE, HCTEM, SLOPEM.THERMM, PERROR , CAL 

61400  2090  FORMAT! IX, /,4X. 'RESPONSE  MAXWELL  1/EX  LETH 

61500  &  THERMAL  PERCENT  CALIB.’,/, 

61600  &  4X, 'MATRIX  TEMP, SHAPE  FACTOR  SLOPE  FACTOR  ERROR 

61700  a  FACTOR’, 

61300  a  /*,  IX , A8 , F9 .2 , ' , '  ,F4 . 2 .F3.4 ,4( F10.4) ) 

61900  2100  WRITE! 1,2110) 

62000  2110  FORMAT! IX/ 10X ,9HDETECT0RS . 10X .8KMEASURED, 5X , 

62100  a  12H  CALCULATED, 5X.7HPERCENT,/ 

62200  X  30X, 6HC0UNTS,  10X ,6HC0t'NTS,6X  ,  10HDIFFERE.NCE) 

62300  DO  2120  1=1, KK 

62400  L=LL( I ) 

62500  2120  WRITE! 1,2130), BALL(L) , BCE! I ) , BCC! I ) , PCTESR! I ) 

62600  2130  FORMAT! 10X , A10, IX ,0PF15. 3.0PF16. 3.0PF15. 3) 

62700  WRITE!  1,2140)  SUMSPC ,  A  VEEN' .  SUMRAD ,  SUMREM ,  QF ,  SUMTLD ,  SLT1HAN , 

62800  X  SUMNTR , SUMNTA , SUMA70 

629 00  2140  FORMAT!/, 

63000  a  10X, 'TOTAL  FLUENCE= ' ,T34, 1PE1 1 , 3.T43 . ' NEUTR0NS/CM2 ' ,/ , 

63100  a  lOX.'AVE  ENERGY(LESS  TH)-' ,T34, 1FE11 .3.T43, 'MEV' ,/ , 

63200  a  10X , ' DOSE= ' ,T34, 1PE1 1 . 3,T43, ' RAD ' ,/ , 

63300  a  10X, 'DOSE  EQUIVALENT-', T34.1PE11. 3, T48, 'REM',/, 

63400  a  10X, 'QUALITY  FACTOR.' .T34.0PF7. 3, T48, 'REM/RAD' ,/, 

63500  a  lOX.'SRLTLD  RESPONSE. ’ ,T34 , F7 . 3, T4S . 'REM/REM(CF-252 )',/ . 

63600  a  10X,  "'HANKINS"TLD  RESPONSE.’ ,T34 ,F7 . 3.T4S, 'REM/REM(CF-252 )’,/ , 

63700  a  10X , 'NEUTRAK  RESPONSE.' ,T34 , F7 . 3, T48, ’REM/REM(CF-252) ',/ , 

63800  a  lOX.'NTA  RESPONSE.' ,T34,F7. 3, T48, 'REM/REM(CF-252) ' ,/, 

63900  a  10X, 'ANPDR-70  RESPONSE. ' ,T34 , F7 . 3 ,T48, ' REM/REM(CF-252 )',/ ) 

64000  WRITE! 1 .2150) 

64100  2150  FORMAT! 3X, 

64200  a  '  BIN  ENERGY  FLLENCE  FLUE.VCE 

64300  a  DOSE  DOSE  EQV.  DOSE  EQV. ' ,/, 

64400  a  '  NO.  MAX  (MEV)  NF.UT/CM2 

64500  a  N7CM2/LETH  (RAD)  (REM)  (C  OF  TOTAL) ' ) 

64600  DO  2160  I-l.JJ 

64700  2160  WRITE! 1,2170)  I ,EEND( 1+1 ) ,SPC( I ) ,SPL( I) ,RAD( I) , 

64800  a  REM(  I )  ,PREM(  I) 

64900  2170  FORMAT! 2X, 14, 2X, 6(1  PEI  1.3)) 

65000  WRITE(1 .2180) 

65100  2180  FORMAT ! 1 X , / , / ) 

65200  2190  CONTINUE 

65300  C  SET  FLAGS,  RETURN  FOR  ANOTHER  SPECTRUM  IF  DESIRED 
65400  JX=0 

65500  LX.O 

65600  CHDET.'N" 

65700  CHNTM.'N" 

65800  2200  WRITE! 5,2210) 

65900  2210  FORMAT!'  LAST  SPECTRUM?') 
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READ(5,1050)LASTSP 
IF(LASTSP.EQ. 'S' )G0  TO  2220 
IF(LASTSP.£Q. *Y’)CO  TO  2240 
GO  TO  2200 

2220  WRITE(5,2230) 

2230  FORMAT!’  CHANGE  MATRIX,  DETECTORS,  UNFOLDING  CODE',/, 
A  '  AND/ OR  ENERGY  INTERVALS? ’ ) 

READ(5 , !050)CHMTX 
IF(CHMTX.£Q. ’N')GO  TO  !220 
IF(CHMTX.EQ. ’Y’)GO  TO  1010 
GO  TO  2220 

2240  IF(FILE.EQ. 'N')GO  TO  2250 

C  WRITE  SPECTRA  TO  FILE  FOR  PLOTTING  IF  REQUIRED 
OPEN ( UNIT- 1 , FILE- ’ SPECX 1 ) 

WRITE( 1,*)SPLPLT 

CLOSE! UNIT. 1 , FILE- ’ SPECX 1 ) 

2250  END 


i 

L 
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APPENDIX  B 

Description  of  Input  Parameters  of  BUNK  I 

SLOPEJ  (9800,38600)  The  initial  value  of  the  slope  of  the  1/E  part  of  the 
MAXIET  spectrum.  Usually  0. 

PERSLP  (9900,37600,28600,38400,41600)  The  amount  by  which  the  slope  of  the 

1/E  part  of  the  MAXIET  spectrum  is  changed  in  searching  for  a  better 
fit  to  the  data.  Typical  values,  0.005  -  0.02. 

THERMJ  (10000,28700)  The  initial  value  of  the  thermal  bin  of  the  initial 
MAXIET  spectrum.  Usually  1.0. 

THMMIN  (10100,41200)  The  minimum  value  allowable  for  the  thermal  bin  of  the 
MAXIET  spectrum. 

THMMAX  (10200,38000)  The  maximum  allowable  value  for  the  thermal  bin  of  the 

MAXIET  spectrum.  TttMMAX  and  THMMIN  are  set  from  physical  characteristics 
of  the  radiation  environment. 

DEAD  (10300,30500)  The  dead  time  of  the  instrument  used  to  determine  the 
detector  counts. 

SHP  (10400,32500,42700)  The  minimum  value  of  the  (1+1)  bin  relative  to  the 
I  bin  for  the  initial  MAXIET  spectrum.  Used  to  limit  the  high  energy 
roll-off  of  the  calculated  Maxwellian  spectrum. 

TSTRAT  (10500,53700)  The  maximum  allowable  value  of  the  error  on  the  fit 
relative  to  the  value  when  the  error  was  last  tested.  Prevents 
further  iterations  when  no  significant  improvement  in  fit  is 
occurring  with  further  iterations.  If  set  >1.0  this  test  will 
not  terminate  the  fit.  Typical  values:  0.9,  0.99,  0.999,  0.9999,  1.1. 

TEMPIJ  (25900-26600,28500,41400)  The  initial  guess  of  the  temperature  of 
the  Maxwellian  peak.  Program  asks  for  "Maxwellian  Temp.". 

SHAPE  (25900-26600,  36400,42600)  The  shape  of  the  high  temperature  portion 
of  the  Maxwellian  peak.  Typical  values,  0-0.5.  Originally  the 
program  searched  for  the  best  shape,  but  this  feature  was  considered 
unnecessary  and  is  now  user  input.  Program  asks  for  "shape". 

PERTMP  (25900-26600,  40100,40200,40500,41400,41500,41800)  The  amount  by 
which  the  Maxwellian  temperature  is  changed  in  searching  for  a 
better  fit  to  the  data.  Should  be  approximately  10%  the  Maxwellian 
temperature,  TEMPIJ.  May  be  positive,  negative,  or  zero.  Positive 
searches  for  a  lower  temperature,  negative  for  a  higher  temperature. 

If  set  to  zero,  the  fit  is  forced  from  TEMPIJ  which  is  often  a 
useful  feature.  Program  asks  for  "perturbation". 

TSTPER  (26900-28000,28400,53800)  An  end  test  used  to  terminate  the  fit. 

When  the  error  on  the  fit  drops  below  this  error  the  fit  is 
terminated.  Program  asks  for  "end  test  (%)". 
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APPENDIX  a  (Con't) 


SMO  (26900-28000,  48800,51900)  The  smoothing  factor.  Typical  values, 
0-0.05.  Smoothing  increases  with  increasing  SMO.  Program  asks 
for  "smoothing  factor". 

CAL  (26900-28000,  56400,60600)  A  calibration  factor  used  to  correct  the 
spectrum  to  agree  with  some  calibration  standard.  Typical  values 
we  have  found  for  a  4  mm  x  4  mm  crystal  and  the  Sanna  matrix  are 
1.2-1. 6  for  Cf-252.  Program  asks  for  "calibration  factor". 

ITRTST  (26900-28000,  48100,51100,53700)  The  number  of  iterations  before 
making  a  test  to  decide  if  the  fit  should  be  terminated.  Typical 
values,  1-100.  Program  asks  for  "iterations  before  error  test". 

ITRMAX  (26900-28000,  45000,47000,53700,54600,59200,59300,60500)  The 

maximum  number  of  iterations  allowed.  Typically  100-1000.  May 
also  be  set  to  0.  If  the  user  inputs  the  initial  spectrum  and 
ITRMAX=0  the  program  does  not  ask  for  detector  data  but  calculates 
the  output  parameters  directly  from  the  input  spectrum.  This 
option  is  useful  for  calculating  integral  parameters  from  known 
spectra.  If  using  the  MAXIET  algorithm,  setting  ITRMAX=0  prevents 
any  further  fitting  of  the  data  by  either  SPUNIT  or  B0N31G.  The 
MAXIET  spectrum,  and  calculated  integral  parameters,  are  then  output 
in  a  slightly  different  format  than  that  shown  in  Appendix  D. 
Program  asks  for  "maximum  number  of  iterations". 

PERTHM  (28100,37900,38500,41700)  The  amount  by  which  the  thermal  bin  of 
the  MAXIET  spectrum  is  changed  in  searching  for  a  better  fit  to 
the  data.  PERTHM  is  not  input  directly,  but  is  calculated  from 
PERSLP . 

PERE  (28200,33000,33200,36700,38400,38500)  The  amount  by  which  the 

magnitude  of  the  1/E  part  of  the  MAXIET  spectrum  is  changed  in 
searching  for  a  better  fit  to  the  data.  PERE  is  not  input  directly, 
but  is  calculated  from  PERSLP.  A  reasonable  relationship  between 
PERSLP,  PERTHM,  and  PERE  must  be  maintained  to  insure  that  the 
MAXIET  algorithm  will  find  the  best  fit  to  the  detector  data. 


28 


APPENDIX  C 

SAMPLE  OF  INTERACTIVE  BUNKI  SESSION 

.LOGIN  720,1564 

JOB  32  RNL-603A-DUAL-KI -54  TTY  127 
[LGNJSP  Other  jobs  same  PPN:63] 

[LGNRDU  Recomputing  disk  usage] 

[LGNQTA  DSKU  IN:3000  0UT.-3000  USED:  1605] 

1510  17-Feb-84  Fri 

[No  mail] 

■EX  BUNKI .FOR 

LINK:  Loading 

[LNKXCT  BUNKI  execution] 

NUMBER  OF  ENERGY  INTERVALS? 

25 

NUMBER  OF  DETECTORS? 

7_ 

TYPE  DETECTOR  CODES  (?  FOR  HELP) 

0_ 

2C 

3C 

5_ 

8_ 

10 

12 

TYPE  MATRIX  NAME  (?  FOR  HELP) 

SAN4 

TYPE  UNFOLDING  CODE  (?  FOR  HELP) 

SPUN  IT 

PUT  SPECTRA  IN  FILE  FOR  PLOTTING? 

N_ 

TYPE  THE  HEADING  (<80  CHARACTERS) 

_ CF  252  in  60  CM  STEEL  BALL  +1.27  CM  LUCITE 

SEARCH  FOR  MAXWELLIAN,  1/E  INITIAL  SPECTRUM? 

N_ 

TYPE  INITIAL  SPECTRUM (26  VALUES) 

26*1 

TYPE:  END  TEST($), 

SMOOTHING  FACTOR, 

CALIBRATION  FACTOR, 

ITERATIONS  3EF0RE  ERROR  TEST, 

AND  MAXIMUM  NUMBER  OF  ITERATIONS 
.2,0, 1,25, 1000 

TYPE  BARE  BONNER  SPHERE  COUNT,  %  ERROR 
935,3 

TYPE  2"+CD  BONNER  SPHERE  COUNT,  %  ERROR 

2125. 1 

TYPE  3"+CD  BONNER  SPHERE  COUNT,  %  ERROR 
5197,  1 

TYPE  5  INCH  BONNER  SPHERE  COUNT,  %  ERROR 

8361.1 

TYPE  8  INCH  BONNER  SPHERE  COUNT,  %  ERROR 

4547. 1 


TYPE  10  INCH  BONNER  SPHERE  COUNT,  %  ERROR 

2305.1 

TYPE  12  INCH  BONNER  SPHERE  COUNT,  %  ERROR 

1063.1 

O.OOOOOOOE+OO,  O.OOOOOOOE+OO,  O.OOOOOOOE+OO, 
O.OOOOOOOE+OO,  0.4940264  ,  1000, 

SAVE  THESE  RESULTS? 

Y_ 

LAST  SPECTRUM? 

N_ 

CHANGE  MATRIX,  DETECTORS,  UNFOLDING  CODE 
AND/OR  ENERGY  INTERVALS? 

N_ 

CHANGE  THE  HEADING? 

N_ 

SEARCH  FOR  MAXWELLIAN,  1/E  INITIAL  SPECTRUM? 

Y_ 

TYPE  MAXWELLIAN  TEMP,  SHAPE,  AND  PERTURBATION 
.5,0, .05 

CHANGE  FIT  PARAMETERS? 

Y_ 

TYPE:  END  TESTC?) 

SMOOTHING  FACTOR, 

CALIBRATION  FACTOR, 

ITERATIONS  BEFORE  ERROR  TEST, 

AND  MAXIMUM  NUMSER  OF  ITERATIONS 
1,0,1,1,1000 
CHANGE  BALL  DATA? 

N 


0.5000000 

,  O.OOOOOOOE+OO, 

0. 1458820 

2.935984 

,  10.29025 

0.4500000 

,  O.OOOOOOOE+OO, 

0. 1255103 

2.985984 

,  8.285000 

0.4000000 

,  O.OOOOOOOE+OO, 

0. 1 118169 

4.299817 

,  5.430347 

0.3500000 

,  O.OOOOOOOE+OO, 

0.8898502E-01 

5.159780 

,  3.348522 

0.3000000 

,  O.OOOOOOOE+OO, 

0.6416974E-01 

6. 191737 

,  1.054216 

0.2500000 

,  O.OOOOOOOE+OO, 

0.9536391 E-02 

2.488320 

,  1.539098 

0.3000000 

,  O.OOOOOOOE+OO, 

0,5167460 

6. 191737 

,  0.7344218 

1, 

SAVE  THESE  RESULTS? 

x_ 

LAST  SPECTRUM? 

Y_ 

END  OF  EXECUTION 

CPU  TIME:  33.33  ELAPSED  TIME:  3:54.08 
EXIT 


O.OOOOOOOE+OO, 


0.1400000 
0.1500000 
0.1800000 
0. 1900000 
0.2100000 
0 . 5000000E-0 1 
0.2100000 
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APPENDIX  D 

SAMPLE  OUTPUTS  OF  BUNKI 


(specification):  b:steve2 

Viewing  file  'b:steve2'.  Hit  <space>  to  continue;  ’Alt-V"  to  terminate. 


TYPE  BUNN . DAT 


CF  252  IN 

60  CM  STEEL  BALL  +  1 

.27  CM  LUCITE 

RESPONSE  UNFOLD 

MAXWELL  CALIB . 

SMOOTH 

PER  CENT  NO.  OF 

MATRIX  CODE 

TEMP, SHAPE  FACTOR 

FACTOR 

ERROR  ITERATIONS 

SAN4  SPUN 

0.00,  .00  1.0000 

0.0000 

0.4940  1000 

DETECTOR 

MEASURED 

CALCULATED 

PERCENT 

COUNTS 

COUNTS 

DIFFERENCE 

BARE 

935.000 

935.163 

0.017 

2"+CD 

2125.000 

2121.371 

-0.171 

3"+CD 

5197.000 

5223.262 

0.505 

5  .INCH 

8361.000 

8307.583 

-0.639 

8  INCH 

4547.000 

4579.329 

0.711 

10  INCH 

2305.000 

2289.550 

-0.670 

12  INCH 

1063.000 

1065.617 

0.246 

TOTAL  FLUENCE- 

3.426E+04 

NEUTR0NS/CM2 

AVE  ENERGY (LESS  TH)= 

3. 104E-01 

MEV 

DOSE= 

4 . 1 15E-05 

RAD 

DOSE  EQUIVALENT* 

3.372E-04 

REM 

QUALITY  FACTOR. 

8.194 

REM/RAD 

NRL  TLD  RESPONSE. 

6.570 

REM/REM(CF-252) 

"HANK1NS"TLD  RESPONSE. 

9.521 

REM/REM(CF-252) 

NEUTRAK  RESPONSE. 

0.521 

REM/REM(CF-252) 

NT A  RESPONSE. 

0.400 

REM/REM(CF~252) 

ANPDR-70  RESPONSE. 

0.992 

RLM/REM(CF-252) 

BIN 

ENERCY 

FLUENCE 

FLUENCE 

DOSE 

DOSE  EQV. 

DOSE  EQV. 

NO. 

MAX  (MEV) 

NEUT/CM2 

N/CM2/LETH 

(RAD) 

(REM) 

{Z  OF  TOTAL) 

1 

4.140E-07 

4.409E+03 

2.727E+03 

2.319E-06 

5.075E-06 

1.5O5E+C0 

2 

6.826E-07 

3.949E+02 

1.818E+03 

2.404E-07 

4.865E-07 

1 .443E-01 

3 

1.445E-06 

3 . 489E+02 

1.071E+03 

2. 154E-07 

4 . 375E-07 

1.29SE-01 

4 

3.059E-06 

2.968E+02 

9.112E+02 

1.821E-07 

3.692E-07 

1.095E-01 

5 

6.476E-06 

2.814E+02 

8,640Ei02 

1 . 708E-07 

3.453E-07 

1.024E-01 

6 

1.371E-05 

2.778E+02 

8. 529E+02 

1.669E-07 

3.364E-07 

9.97SE-C2 

7 

2 . 902E-05 

2.939E+02 

9.024E+02 

1.755E-07 

3 . 512E-07 

1.042E-01 

8 

6. 144E-05 

3 . 167E+02 

9.722E+02 

1.880E-07 

3 . 734E-07 

1 . 107E-01 

9 

1.301E-04 

4.070E+02 

1.249E+03 

2. 398E-07 

4.726E-07 

1 . 402E-01 

10 

2.754E-04 

3 .  S89F.+02 

1 . 194E+03 

2.216E-07 

4 . 368E-07 

1.295E-01 

11 

5.929E-04 

4 . 741E+02 

1.424E+03 

2.591E-07 

5. 121 E— 07 

1.51«E-01 

12 

1.234E-03 

5. 177F.+02 

1.626E+03 

2. 718E-07 

5.384E-07 

1 . 597E-01 

13 

2.613E-03 

6. 231E+02 

1 . 912E+03 

3.208E-07 

6.356E-07 

1.8S5F.-01 

14 

5.531E-03 

7 . 779E+02 

2.3S9E+03 

3.954E-07 

7.S42E-07 

2.32oE-01 

15 

1 . 171E-02 

9.992E+02 

3.067E+03 

5.035E-07 

1 .00SE-06 

2.990E-01 

16 

2.479E-02 

1 . 324E+03 

4.066E+03 

7 . 456E-07 

2.005E-06 

5.947E-01 

17 

5.247E-02 

1 .S70E+03 

5.741E+03 

1 .241E-06 

5.023E-06 

1 . 490E+00 

13 

1.11  IE-01 

2.903E+03 

8.909E+03 

2.27SF.-06 

1 . 386E-05 

4.111  F.+uO 

19 

2.237F.-01 

4.483E+G3 

1 .475E+04 

4.654E-06 

3.721E-05 

1 . 104F.+01 

20 

4 . 508E-0 1 

6.445E+03 

2.118F.+04 

9.371E-06 

9. 126E-05 

2.7O7E+01 

21 

9.072E-01 

4.017E+03 

1 . 520E+04 

1.036E-05 

1 . 104E-04 

3.274E+01 

22 

I -872E+00 

1.422E+03 

4. 520E+03 

4.844E-06 

4 . 996E-05 

1 . 479E+01 

23 

3.679E+00 

2 .8S6F.+02 

9.834E+02 

1 . 218F.-Ob 

1 . 1 53E-05 

3.419F.+G0 

24 

7.408E+00 

7.049E+01 

2.319E+02 

4.073E-07 

2.871 E-06 

8.513E-01 

25 

1.492E+01 

2.413E+01 

7.935E+01 

1 . 596E-07 

9.919E-07 

2.942E-01 

X 

A 
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1  CF  252  IN  oO  CM  STEEL  BALL  +  1.27  CM  LL'CITE 


RESPONSE  UNFOLD 

.MAXWELL  CALIB . 

SMOOTH 

PER  CENT  NO.  OF 

MATRIX  CODE 

TEMP. SHAPE  FACTOR 

FACTOR 

ERROR  ITERATION’S 

Sana  spun 

0.30,  .00  1.0000 

0.0000 

0.7344  1 

DETECTOR 

MEASURED 

CALCULATED  PERCENT 

COUNTS 

COUNTS 

DIFFERENCE 

BARE 

935.000 

931.928 

-0.329 

2"+CD 

•  2125.000 

2105.156 

-0.934 

3"+tD 

5197.000 

5234.559 

0.723 

5  INCH 

8361.000 

8283.147 

-0.931 

8  INCH 

4547.000 

4596.901 

1.097 

10  INCH 

2305.000 

2303.379 

-0.070 

12  INCH 

1063.000 

1067.717 

0.444 

TOTAL  FLl'ESCE. 

3.420E+04 

NEUTRCNS/CM2 

AVE  ENERCY( LESS  TH)» 

2.94SE-01 

MEV 

DOSE- 

4 . 106E-05 

RAD 

DOSE  EQUIVALENT* 

3.396E-04 

REM 

QUALITY  FACTOR* 

8.270 

REM/RAD 

NRL  TLD  RESPONSE* 

6.597 

REM/REM(CF-252) 

"HANKINS"TLD  RESPONSE* 

8.929 

REM/REM(CF-252) 

NEUTRAK  RESPONSE* 

0.492 

REM/REM(CF-252) 

NTA  RESPONSE* 

0.399 

REM/REM(CF-252) 

A.NPDR-70  RESPONSE* 

0.999 

REM/REM(CF-252) 

BIN 

ENERGY 

FLUENCE 

FLUENCE 

DOSE 

DOSE  EQV. 

DOSE  EQV. 

NO. 

MAX  (MEV) 

NEUT/CM2 

N/CM2/LETH 

(RAD) 

(REM) 

(Z  OF  TOTAL) 

1 

4. 140E-07 

4 . 709E+03 

2.912E+03 

2.477E-06 

5.420E-06 

1 . 596E+00 

2 

6.826E-07 

1.015E+02 

4.675E+02 

6 . 180E-08 

1.251E-07 

3.683E-02 

3 

1.445E-06 

1.732E+02 

5.317E+02 

1.069E-07 

2 . 172E-07 

6.394E-02 

4 

3.059E-06 

2. 024 £+02 

6.215E+02 

1 . 242E-07 

2.513E-07 

7.415E-02 

5 

6.476E-06 

2.367E+02 

7.267E+02 

1.437E-07 

2.904E-07 

8. 552E-02 

6 

1.371E-05 

2 . 768E+02 

8.498E+02 

1.563E-07 

3.352E-07 

9.87] E-02 

7 

2.902E-05 

3.236E+02 

9.938E+02 

1.932E-07 

3.868E-07 

1.139E-01 

8 

6 . 144E-05 

3.787E+02 

1.162E+03 

2.248E-07 

4.464E-07 

1 .315E-01 

9 

1.301E-04 

4. 430E+02 

1 .360E+03 

2.610E-07 

5. 144E-07 

I . 515E-01 

10 

2.754E-04 

5.181E+02 

1.591E+03 

2.952E-07 

5.819E-07 

1.713E-01 

11 

5.929E-04 

6.211E+02 

1.S65E+03 

3.394E-07 

6 . 708E-07 

1.975E-01 

12 

1.234E-03 

6.953E+02 

2.184E+03 

3.651E-07 

7 . 232E-07 

2. 129E-01 

13 

2.613E-03 

8. 335E+02 

2 . 55SE+03 

4.292E-07 

8.501 E-07 

2.503E-01 

14 

5.531E-03 

9.817E+02 

3.014E+03 

4.990E-07 

9.S93E-07 

2.914E-01 

15 

1 . 171 E— 02 

1 . I 70E+03 

3. 592E+03 

5 . 895E-07 

1.180E-06 

3.476E-01 

16 

2.479E-02 

1.431E+03 

4.393E+03 

8.055E-07 

2 . 167E-06 

6.380E-01 

17 

5.247E-02 

1.850E+03 

5.630E+03 

I.22SE-06 

4 . 970E-C6 

1.463E+O0 

13 

1.111E-01 

2.615E+03 

8.025E+03 

2.052E-06 

1 .249E-G5 

3.676E+00 

19 

2.237E-01 

3.653E+03 

1 . 202E+04 

3.792E-06 

3.032E-05 

8.927E+00 

20 

4 . 508E-01 

5.633E+03 

1.S51E+04 

8. 191E-06 

7.977E-05 

2.349E+01 

21 

9.072E-01 

5.486E+03 

1 . 806E+04 

1 . 231E-05 
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